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We study the effects of Coulomb interaction between 2D Weyl fermions with anisotropic disper¬ 
sion which displays relativistic dynamics along one direction and Newtonian dynamics along the 
other. Such a dispersion can be realized in phosphorene under electric field or strain, in Ti02/V02 
superlattices, and, more generally, at the quantum critical point between a nodal semimetal and an 
insulator in systems with a chiral symmetry. Using the one-loop renormalization group approach in 
combination with the large-A expansion, we find that the system displays interaction-driven non- 
Fermi liquid behavior in a wide range of intermediate frequencies and marginal Fermi liquid behavior 
at the smallest frequencies. In the non-Fermi liquid regime, the quasiparticle residue Z at energy 
E scales as Z oc with a > 0, and the parameters of the fermionic dispersion acquire anomalous 
dimensions. In the marginal Fermi-liquid regime, ^ oc (| logi?|)“*’ with universal b = 3/2. 


Introduction .—After the discovery of time-reversal in¬ 
variant topological band insulators mm, the notion of 
topological states of matters has been extended to a 
broad class of systems. In particular, recent studies of 
three-dimensional (3D) Weyl and Dirac semimetals with 
nodal points have demonstrated that these systems also 
possess quantized topological invariants (or topological 
charges) and associated topological surface states [3HTT]. 
Since the topological invariant assigned to each nodal 
point guarantees its stability, the transition from a topo¬ 
logical semimetal to an insulator can be achieved when 
pairs of nodal points with opposite topological charges 
merge at the same momentum |12] . Hence the quantum 
critical point (QCP) of semimetal-insulator transitions 
should have emerging gapless degrees of freedom with 
zero topological charge umii. 

Because the topological charge of a nodal point is solely 
determined by the energy dispersion around it, the van¬ 
ishing of a topological charge at a QCP implies that at 
this point the dispersion of low energy excitations must 
become unconventional. Indeed, it has recently been 
shown that a new type of fermionic excitations, dubbed a 
3D anisotropic Weyl fermion (AWF), appears at the QCP 
between a 3D Weyl semimetal and an insulator 
A 3D AWF has an anisotropic energy dispersion, which 
is quadratic in one direction and linear in the other two 
orthogonal directions. Such a dispersion brings about 
highly unusual quantum critical behavior. Most notably, 
quantum fluctuations of 3D AWFs screen the long-range 
Coulomb interaction and make it anisotropic, however 
the long-ranged nature of the interaction is preserved. 
The screened anisotropic Coulomb potential becomes an 
irrelevant perturbation in the low-energy limit, i.e., low- 
energy fermions remain free quasiparticles Ca¬ 
in this letter, we describe a semimetal-insulator tran¬ 
sition and associated quantum criticality in a system of 
two-dimensional (2D) AWF with long-range Coulomb in¬ 


teraction. The role of the Coulomb interaction in 2D 
nodal semimetals V{q) oc l/|q| has been widely studied 
with particular emphasis on graphene |15H18j . The con¬ 
clusion is that at strong coupling, the Coulomb interac¬ 
tion generates anomalous exponents |18] . This behavior, 
however, holds only at intermediate frequencies because 
the dimensionless coupling flows towards smaller values, 
and below a certain energy the system necessary enters 
into a weak coupling regime. In this regime the renor¬ 
malizations are only logarithmical (marginal) [TfllTM^ . 
and the quasiparticle residue tends to a finite value at 
zero energy, i.e., at smallest frequencies graphene pre¬ 
serves Fermi-liquid behavior. As a consequence, interac¬ 
tions dress physical observables, like the optical conduc¬ 
tivity, only by extra logarithmic factors [HlEolEllEe]. 

We argue that the behavior changes fundamentally 
when a semi-metal is brought to the quantum criti¬ 
cal point of the semimetal-insulator transition, at which 
pairs of nodal points merge. In this case, the low en¬ 
ergy excitations around a gapless point are 2D AWF with 
linear dispersion in one direction and quadratic in the 
other [27li35] . i.e., a 2D AFW displays simultaneously 
relativistic and Newtonian dynamics. 

We analyze the effects of Coulomb interaction in 2D 
AWF’s by combining a renormalization group (RG) anal¬ 
ysis and a large N expansion. We find (see Fig. that 
over a wide range of energies the screened Coulomb in¬ 
teraction is a relevant perturbation, and the system is 
in the strong coupling limit and displays non-Fermi liq¬ 
uid behavior with power-law energy dependence of the 
quasiparticle residue Z oc E°', where a = 0{{}ogN)/N). 
This behavior starts at Ei ~ '^ogN^ where A is of 

the order of the bandwidth and bi = 0 ( 1 ), and extends 
down to very low energy E^ ~ A/{N^). In the subrange 
E 3 < E < E 2 < El, where E 2 ^ 62 = 0(1), the 

parameters of the fermionic dispersion (the velocity and 
the effective mass) also become energy dependent and 
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FIG. 1: (Color online) Evolution of quasiparticle properties 
as the energy scale varies. There are three energy scales Ei ~ 
^^-b-iN/io^N^ i ?2 ~ , and E 3 ~ A/N^ characterizing 

each region. Here foi ,2 are constants of 0(1), and v{E), A{E), 
Z{E) are the velocity, the inverse mass, and the quasiparticle 
residue, which are running as E decreases. A is of the order 
of the band width. 



FIG. 2: (a) Feynman diagrams representing the RPA boson 
propagator (bold wavy line). Each fermion loop, accompanied 
by a factor N, is resummed in the boson propagator, (b) 
electron self-energy and (c) the vertex function. 


vary as powers of E with anomalous exponents 0{1/N). 
At even smaller energies E < E 3 , the system crosses 
over to weak coupling behavior. However, contrary to 
the case of 2D nodal semimetals, 2D AWF do not be¬ 
come free particles at the smallest frequencies. Instead, 
in the limit A —f 0, the system displays marginal Fermi 
liquid behavior with universal, Af-independent quasipar¬ 
ticle weight Z{E) oc (log 

The model. —Non-interacting 2D AWF are described 
by 

Ho = -AOItx - ivdyTy, ( 1 ) 

where 1/(2A) > 0 is the mass along x direction and v is 
the velocity along y direction. The Pauli matrices Tx^y 
are used to denote the valence and conduction bands. 
The absence of Tz in the Hamiltonian ensures the chiral 
symmetry of the system, thus a nodal point can have an 
integer winding number. The semimetal-insulator tran¬ 
sition can be described by adding to Hq a perturbation 
term rriTx with a constant m. Depending on the sign of 
m, the system becomes either a gapped insulator (m > 0) 
or a semimetal {m < 0) with two nodal points on the kx 
axis with the winding numbers ±1, respectively. The 
distance between two nodal points decreases as |m| is re¬ 
duced, and at m = 0 two nodal points merge, resulting 
in a gapless point with a zero topological charge [371HH]- 
The anisotropic dispersion leads to the density of states 
PQGp{E) oc -y/A, which is obviously enhanced in the low 
energy limit as compared to the case of a semimetal for 
which psm{E) oc E. This suggests that interaction effects 
are most relevant at the QCP. 

We consider long-range Coulomb interaction between 
2D AWF. The corresponding effective action is 

S= J drd'^xil^lKdr+ig(l)) + Ho]^|^a + ^ J drd^x{di'p)‘^, 

( 2 ) 


the summation over repeated indices is implied. Ob¬ 
serve that the bosonic field (p is defined in 3D space 
whereas the electron is confined to a 2D plane. Once 
the z-dependence of (j>, is integrated out, the Coulomb 
potential in momentum space becomes V{q) oc l/|q|. 
The dimensionful boson-fermion coupling associated with 
Coulomb potential is g = ej^fe., where e is the electric 
charge and e is determined by the dielectric constant. 
The corresponding dimensionless coupling a, which ap¬ 
pears in perturbation theory, is the ratio of the Coulomb 
potential E^ ^ A~^vg‘^ and the electron kinetic energy 
Ekin A~^v‘^: a = = ^- To control the theory 

analytically, we extend tEe model to N fermionic flavors 
and consider the large N limit. At large N, the dimen¬ 
sionless coupling constant becomes oat = Na. 

Bosonic and fermionic propagators .—We follow the 
same strategy as in earlier approaches on large N theories 
of quantum-critical behavior of itinerant fermions mm 
and compute fermionic and bosonic self-eneries in a self- 
constent scheme. Namely, we first compute the one- 
loop bosonic self-energy (bosonic polarization operator 
ni(D, q)) as the latter contains N (and also, as we will 
see, contains smaller power of Qy compared to the bare 
term then use the dressed dynamical 

bosonic propagator 

Df\n,q)^D^\q)-NU,{n,q), (3) 

to compute the fermionic self-energy and corrections to 
the one-loop ni(q, D) within the 1/A^ expansion. We 
show that the corrections depend logarithmically on the 
running energy and solve the RG equations for the full 
propagators using 1/A^ as a control parameter. 

The evaluation of the one-loop polarization bubble is 
rather involved. We present the details in the Supple¬ 
mental Material (SM) and here list the result. We find 
that ni(0, q) can be expressed as 


where ipa describes a two-component fermion field with 
the subscript a = 1,...,N labeling the species of 
AWF, and ^ is a bosonic field which one obtains via 
Hubbard-Stratonovich transformation of the instanta¬ 
neous Coulomb potential. The subscript i = x,y,z and 


ni(H,q) = -a 


dxA^I'^gl dyA 
A(D,q)i/4+ A(0,q)3/4 ’ 


( 4 ) 


where A(r2, q) = and dx, dy, and c are 

constants, whose explicit values we present in SM. This 
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Hi {il, q) matches the exact results, which we obtained 
analytically, in the three limits: (i) Hi (17 = qy = 0 ) = 


(ii) ni(|n| > Aql,qy = 0 ) = 


g r(3/4) VAql 
8V^r(9/4) yn ’ 


(hi) ni(g, = 0) = Observe that 

while Hi(17 = qy = 0) oc \qx\ has the same functional 
dependence as DQ^{qx) = \qx\, T{i{qx = 0,17 vqy) oc 
is parametrically larger than DQ^{qy) = \qy\ at 
small qy (Ref. |37j). 

We now use the bosonic propagator with ni(17,q) in¬ 
cluded and compute the one-loop fermionic self-energy 
Ei(a;, fc) and vertex correction 5gi. The corresponding 
diagrams are shown in Fig. For the self-energy, we 
obtain at large oat = Ng'^/v (see SM for details.) 


Ei{u},k) 


= {-i9f j ^|^|Go(a;-|-17,fc-f q)I7(17,q) 

= '^UJ ‘ ‘ Ak^Tx '^ky ' '^kyTyj (5) 


in which 


= Izl, + 'JaI, ^ky = + Ivl, ( 6 ) 

where I = log{A/E), A is the upper energy cutoff, 
and E is the largest of {\uj\,VF\ky\, Ak"^). The param¬ 
eters are log a at/A^, = 0.1261/A^, and 

y^ = 0.3625/iV. We see that the quasiparticle residue 
Z = (1 + dE/diuj)~^ and the parameters of electronic 
dispersion acquire logarithmically singular 1/iV correc¬ 
tions. Observe that the correction to Z is stronger than 
the corrections to v and A by logaAr. 

For the vertex correction we find at vanishing external 
momentum and frequency 

/ d^k 

— Dik)Goik)Goik) = Eyy. (7) 

This is consistent with the Ward-Takahashi identity. The 
calculation of Sgi generally requires care as Sgi has terms 
which depend on how the limit of zero momentum and 
frequency is taken |38j . However these terms are non- 
logarithmical and can be safely neglected at large N. 

Renormalization group analysis. —One can verify that 
higher-order corrections contain higher powers of h To 
sum up the series of logarithms, we can either use the 
Wilsonian shell RG analysis [3S], or express the full self¬ 
energy S(a;, k), the full vertex g, and the full polarization 
bubble self-consistently, via full Green’s functions and 
full vertices, and then represent E{1), g{l), and n(^) as 
integrals f dl' over running I', and obtain RG equations 
by taking derivatives with respect to the upper limit m- 
Either way we obtain the same set of RG equations for 
running Z, v, and A. We also verified that the product 
gZ is not renormalized, as it is required by the condition 
that the electric charge is a conserved quantity. The RG 
equations for Z{1), v{l), and A{1) are {X = dX/dl) 


Z{1) = -lz{l)Z{l), v(l) = y„u(0, A(l) = jaAH) (8) 


where yz(0 = yz + y^ log Solving these equations 
and using I = \og{A/E) we obtain v{E), A{E), and Z{E) 
at energy E in the form 


V 


Ay- A{E) 

e) 




-73 + 


-y/TS 'Yv 
^ 3/2 N 


I 


(9) 


The analysis of Eq. ([^ shows that there are three energy 
scales characterizing the system’s behavior. At high en¬ 
ergies, E > El = where bi = 0(1), the 

dependence of the fermionic propagator on E is weak, 
i.e., fermions behave as almost free quasiparticles. At 
E 2 < E < El, where E 2 = Ae~^^^, 62 = 0(1), v 
and A remain close to their bare values, but the quasi¬ 
particle residue becomes strongly E-dependent, and the 
fermionic propagator at the typical energy E acquires a 
non-Fermi-liquid form with anomalous exponent y^, i.e., 
G oc 1/E^“^*. This behavior holds also at energies below 
E 2 , but now V and A grow as powers of A/E with anoma¬ 
lous exponents y„ and y^i, respectively. The presence of 
anomalous dimensions in the theory implies that physi¬ 
cal observables, such as the specific heat, the compress¬ 
ibility, the diamagnetic susceptibility, etc., show unusual 
temperature dependencies, as shown below. 

The strong coupling results differ quantitatively but 
not qualitatively from the case of graphene. In both 
cases. Coulomb interaction gives rise to anomalous ex¬ 
ponents for the quasiparticle residue and the fermionic 
dispersion. 

Weak coupling limit .—The behavior described by 
Eq. (§ holds as long as the dressed dimensionless cou¬ 
pling otN = {gZ)'^N/v remains large. The bare value 
of OAf is of order N, however v{E) grows upon the flow 
towards lower energies, and eventually, at E < E 3 = 
A/N^, the dimensionless coupling a at becomes smaller 
than one. Once this happens, the RG equations have to 
be modified because, e.g., bare \qx\ in the bosonic prop¬ 
agator becomes larger than the \qx\ term in the polar¬ 
ization operator. We evaluated the one-loop self-energy 
and vertex corrections at small un and found that they 
are again logarithmical, but the RG parameter I is now 
I = log(Em/E), where Em = g'^N'^/A is the energy scale 
below which the I term in the polarization domi¬ 
nates over the bare term |gj,| in the boson propagator. 
The prefactors y^, y„ and y^ are different from those in 
Eq. ([^ and are given by 


_ 3 aN _ _ cxn __ _ aw I log Oat I 

“ STT^iV’ ~ 2n^N ’ 


( 10 ) 


where oat = g'^N/v and gZ is not renormalized. Higher- 
order corrections give higher powers of logarithms, and 
performing the same analysis as in the strong coupling 
limit, we hnd that Z{1), v{l), and A{1) satisfy the same 
RG equations as in the strong-coupling limit, but with 
y’s from Eq. (10). Solving these equations we obtain at 
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FIG. 3: (Color online) RG flow of the coupling constant a oc 
1/w, the effective mass 1/A, and the quasiparticle residue Z 
for 2D AWF. We set the initial values aw = 4 and A = \ 
and the number of flavors N — A:, a, 1/A, Z all flow to zero. 
The quasiparticle residue has power-law dependence Z{E) oc 
E°' at intermediate energies, typical for a non-Fermi liquid, 
and scales as Z{E) oc (1/| log at the lowest energies, 

i.e., at vanishing E the system displays marginal Fermi liquid 
behavior. 


smallest energies (largest 1 ) 

2 

vil) = -^l, Z{1) = Z-3/2, A(0 = Aei° 8 '(11) 
47r^ 

We see that the quasiparticle Z does not reduce to a 
constant in the limit Ojv —t 0 but keeps decreasing, even 
at the smallest E. At the same time, we see that the 
series of logarithms at weak coupling do not generate 
anomalous dimensions, i.e., v{l) and Z{1) behave as pow¬ 
ers of I (the inverse mass A has a somewhat more com¬ 
plex dependence on I, but still there is no anomalous 
dimension.) The logarithmic form of Z{1) implies that 
the fermionic Green’s function at the smallest energies 
behaves as G{E) ^ {logE)~^^‘^/E. Such behavior is gen¬ 
erally termed as a marginal Fermi liquid. 

This last result shows that the weak coupling behavior 
of AWF is qualitatively different from that in graphene. 
In graphene, the quasiparticle Z factor tends to a finite 
value at zero energy m, i-e. the system retains Fermi- 
liquid behavior with well defined quasiparticles. AWF, on 
the contrary, do not become sharp quasiparticles, even at 
the lowest energies. 

To verify our analytical analysis, we obtained the RG 
flow of Z, Oat oc 1/v, and 1/A numerically for = 4. 
(See Fig. §) We clearly see that v and A increase upon 
the system flows to lower energies (higher 1 ), whereas the 
quasiparticle residue Z decreases, initially by a power 
law, and then nearly flattens at the largest 1 . 

Physical observables.- Now we obtain scaling relations 
for physical observables. In general, an operator O with 
the scaling dimension zq obeys 

O (k, a;, V, A) = zqO (k (1) ,oj {1) ,v (1) ,A (1)) , (12) 

where k(Z) = {zx{l)kx,e^ky) and Lull) = z^{l)oj are 
the running momenta and frequencies. Zx,ui are deter¬ 
mined by = 1 - 7 A and = 1 + 7 ^ - 7 „. In 


the case of the particle density O = n, zo=n = 
Standard scaling arguments then yield, at strong cou¬ 
pling, for the compressibility k = oc 

4> = + \lA ~ 0.4255/A^. Similarly, the heat ca¬ 
pacity C oc and the diamagnetic susceptibility 

Xdia CX -T- 2 - HU. Particularly interesting behav¬ 
ior is expected for the anisotropic transport and opti¬ 
cal properties. The real part Ua (w) of the optical con¬ 
ductivity for non-interacting AWF is anisotropic with 
'^x,y (^) (w/wo)^^^^, where wq = w^/A and the up¬ 

per (lower) sign stands for the x (y) direction. Using 
gauge invariance, one finds that the scaling dimensions 
of the optical conductivity are different for the two direc¬ 
tions: Za^ = ZxC~^ and = z~\ In the strong coupling 
regime, we immediately obtain that 

o-x.y (w) oc f — j (13) 


with (j)^ = — 1a/‘2- ~ 0.299/A^. Thus, the anisotropy 

of the optical conductivity is amplified by strong inter¬ 
actions. In the weak coupling limit holds instead: 


ax,y (w) oc N— 


e- I a;e'°s i 




Wo log^ 5 


(14) 


Conclusion .—In this paper we studied quantum criti¬ 
cal behavior of 2D interacting AWF. We showed that the 
interplay between the anisotropic electron dispersion and 
the long-range Goulomb interaction generates a highly 
anisotropic structure of the screened Goulomb potential, 
which, in turn, induces unconventional behavior of elec¬ 
trons. We demonstrated that interacting 2D AWFs dis¬ 
play non-Fermi liquid at intermediate energies, with var¬ 
ious anomalous physical properties, and marginal Fermi 
liquid behavior at the smallest energies. 

There are several candidate materials for 2D AWF. 
In deformed graphene , pressured organic conduc¬ 

tor a-(BEDT-TTF) 2 l 3 [H], and artificial lattices of cold 
atoms |42] , 2D AWF emerge via the merging of two Dirac 
points. In Ti02/V02 nanostructures 2D AWF were pre¬ 
dicted to be intrinsic low energy excitations due to the 
peculiar symmetry of the system [5TU55] . 2D AWF were 
also predicted theoretically to exist under external elec¬ 
tric field or strain in black phosphorous (a system with a 
few layers of phosphorene) [301 US] • This prediction has 
been confirmed in a recent angle-resolved photo-emission 
study |44j . We believe that intriguing non-Fermi liquid 
physics can be probed in these systems through a sys¬ 
tematic investigation of quasiparticle properties. 
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Note added .—After the completion of our paper, we 


became aware of the independent work by Gil Young Cho 
and Eun-Gook Moon on the same problem. 
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Supplemental Material 

The effective action describing the two-dimensional anisotropic Weyl semimetals in momentum space with Matsub- 
ara frequency is 


where 


Sf = 


duj (Pk 

(27r)^ 


S = SF+SB+Sg, 


ijpk){—iu} + Ak^Tx -b vkyTy)ilj{u}, k), 






J 2tt 27r (27r)^ (27r)^ 

The fermion propagator Go{ijJ, k) and the boson propagator reduced to the 2D plane Do{q) are given by 

Go(w, fe) = 


Do{q) = J 


— iiO + Ak'^Tx + VkyTy ’ 


dk. 


1 


1 


27r kx + k'^ + k1 




The vertex of fermion and boson propagators gives {—ig). 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 


I. POLARIZATION FUNCTION 


The polarization at one-loop level is given by 


U{n,q) = J ^ J ■^^^Ti[Go{n + u},q + k)Goiuj,k)] 

duj f d?k —uj{Vt + u)) + A^k^^kxqxY'+ v'^ky{ky + qy) 


= tt- 


27r J (27r)2 [(U -|- w)^ -|- A'^{qx + kxP + v‘^{qy -b kyP][uj‘^ -b A^k'^ + v'^k'^]' 

By using Feynman parameterization 

J_= ^_ 

AB Jo [xA+il-x)Bp^ 

we can perform the w and ky integrations analytically to obtain 


1 = 


9^ P , -“^xil - x)qy + kl{kx + qx)^ - x{kx + qx)^ - - x)kx 

-= / dkx / dx - - - 

iir'^vy/A J-oo Jo 


a;(l - x)n'^ + x{l - x)ql + x{kx + qxP + (1 - x)k‘^ 


( 21 ) 


( 22 ) 


(23) 
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A. = 0 and = 0 


Assuming = qy = 0 we have 


n(r! = o,^,^ = iogO) = 


o2 /*oo /•! 7,2 


dkrr, / dx 


47r‘^v^/A J-oo Jo 


kl{kx + - x{kx + - (1 - 

x{K + QxY + (1 - x)k% 


J-oo " I [{kx + qxY - kt\ [ kt 


Using the relation 


j X‘^{X + 1)^ , [(XA 1)4] ) _ TT^ 

[(cc + l)4-a;4] ^4 Jj 4> 


we obtain 


II(f2 — 0,qx, Qy — 0) — 16 ^' 


B. B ^ Aq^ and Qy = 0 


For U ^ Aq"^ and qy = 0, we consider the lowest order contribution in (Ag^/fl), and then the polarization becomes 

Tj {O 9^ 9y _ _ 9^ ,, , kl{k^ + - x{k^ + q^)^ - (1 - x)kt 

V^VA^v J Att^vVaJ-oo Jo xil-x)n^+xik^ + q,r + il-x)kt 
9^ ql [^ , l-12x0-x) 

4:y/2TTv \/M} J0 ^ [a;(l-x)]4/4 


Here B{a, b) is the beta function and therefore we obtain 


OL r(3/4) '/A(^^ r 

n(U,g,,g,- 0 )-- 8 ^r( 9 / 4 ) ^ ^ 


C. (Ji, = 0 

For q^ = 0, we can calculate the polarization function including the frequency dependence: 

= 0,^) = [ dk, f\x Ml 

V \/A vj dir'^vy/Aj Jo [xO - x){n^ + q^) + k*] 


f dxx4/4(l-a;)i/4, 
Jo 


2V27TvVA{n^+q^yr/^ Jo 

9^ 9l p A 5\ 

(U2 + <72)3/4 \a'a)- 


Therefore the polarization for < 7 ^: = 0 becomes 


n(U,<7.=0,<7,) = -^^^ 


a r(5/4) v'^qllVA 
G-v/tt r(3/4) (U 2 _|_ 




FIG. 4: Comparison of the ansatz Fi{u) (blue line) with the exact function ^ 2 ( 11 ) (red dots), which is numerically evaluated 
at each point. The horizontal axis measures u = 4:fl/{Aq^). 


D. An ansatz for the polarization fnnction 


We propose the following ansatz for the polarization function: 

dxA^/'^ql 


n(f^i qx, qy) — 


dyA 


(n^ + ri2g2 _j_ cA 2(74)1/4 (f]2 _|_ y2q2 _|_ cAP-q^'iA ’ 


where 


dx — 


1 r(3/4) 1 r(5/4) _ / 2 r(3/4) y 

8v^r(9/4)’ 8v/)fr(7/4)’ 


(31) 


(32) 


This ansatz is consistent with the exact result of 11 in the following three limits; i.e., (i) 14 = 0 and qy = 0, (ii) 
14 Aq^ and qy = 0, and (iii) qx = 0. 

Since 11(14, q^ = 0, qy) is consistent with the exact result, let us check the behavior of 11(14, qx,qy = 0) more carefully 
and compare it with the exact result: 


where 


fl{n,qx,qy = 0) = alga;I Fj 


1 .rW 4 ) 




Sy^TT r(9/4) {u + c)4/4 

On the other hand, the exact formula of the polarization for = 0 is 

/ 414 

n(14,ga,,gy = 0) = -a|ga:|J^2 f ^ 

where 


_L f rlxrl (l/(^+ 1)"^ + a/(^- i-)^ + p'^)[iz^ - 1)^ {z + 1)^ + p'^ {z - 1)^ + 

StT^ J ^ \/iz + I)"* +P^\/{z - 1)^ +p^[v? + (-\/{z + 1)^ +p^ + \J{z- 1)4 +p 2 ) 2 j 


V(^-l)4+p2)2] 


(33) 


(34) 


(35) 


(36) 


Fig. ^shows Fi{u) (ansatz) and F 2 {u) (exact) as functions of u = 414/(Ag^). We confirm that the ansatz Fi{u) is in 
good agreement with the numerically-estimated exact function F 2 {u). 


II. SELF-ENERGY AT LARGE un 


By using the ansatz polarization function, the inverse boson propagator becomes 


D ^(414, q) = Z4g ^(114, q) — 7Vn(114, q) 


= 2 


\/qx+qy+0:N 


dxA^I\l 


dyA 4/2u2q,2 


(142 


u2g2 


cAFq^^A (1]2 _|_ _|_ g^2g4)3/4 


(37) 
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where = Na, and the electron self-energy is 


1 f d'^ 

T,{iuj,k) = {-ig)‘^-Y] / -^Go{iuj + in,k + q)D{in,q) 

(Pq i{uj -I- fl) -I- A{kx + qx)‘^Tx + v{ky + qy)Ty 




in 

\ _ r rl‘^^ 


^ m {uj+ n)^ + A‘^{kx + qxY+ v^{ky + qyY 


D{in,q), 


(38) 


The self-energy has a constant part S(a; = 0, fe = 0) « —a/{2Tr'^)vA\og{AA/v)Tx, It does not have a divergence in 
the infrared limit, and hence this correction can be neglected in the low-energy analysis. Also it does not alter the 
discussion in the low-energy region. At T = 0, the sum of Matsubara frequency turns to be integral of frequency 
along the imaginary axis. We introduce and to simplify the presentation as follows: 


S(a;, fe) = y ^ J ■^J^Goiiuj + in,k + q)D{in,q) 

= Y.^-iui - • Ak^Tx - ^ky ■ vkyTy. 


(39) 


A. Evaulation of self-energy at T = 0 

1. Soj 


Under kx = ky = 0, we assume |n|, Eq |w| where Eq = ^A^q^ + v'^qy, and expand Go(w -I- O, q) with respect to 
u! up to the linear order. Then we obtain 




dfi r (Pq n^-Ei 

^ J (2^)2 (r!2 + Ely 


D{n,q). 


(40) 


If D{V,,q) is independent of ft, e.g., if we adopt Do{q), the integral over H yields zero because the integrand has no 
first-order pole. We introduce dimensionless variables x and y as 


VAqx _ vqy 


to rewrite the integral. Then we have 

E - 
" (2^)3 

where E^°^ and eD^ are defined by 

e(°) = 

and the function d{Vl,x,y) is given by 




1 — 

(1 + + y‘^Y 


a 

r dD 

- (2^)3 J 

' M 

a 

r dfi 

~ (2^)3 j 

y 

' M 


d(fI,x,y) = E(o)+E«, 


d{il,x,y), 


(41) 


(42) 


(1 -I- -I- 2/2)2 

—x'^ — 2/^ 
^(1 -I- -I- 2/^)2 


di^,x,y), 


d (n, X, y) = 2\ x"^ + 


\n\Ay^ 


■ OiN 


dxx"^ __ dyy'^ 

{1cx'^ E y'^y/'^ {lE cx'^ + y'^Y/‘^\ 


(43) 


(44) 


Note that d{El, x, y) is singular at x = y = Q even for un 1. Due to this singularity, E^^^ and E^)^ will be calculated 
separately. 

To calculate E^^^, we introduce the polar coordinate x = r cos 9, y = r sin 9. In this case, the first term of d ^(17,0;, y) 
must be retained to avoid the singular behavior at r = 0. Dominant contribution comes from small r region. For D 
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integral, we introduce the upper energy cutoff A and the lower cutoff E. Then we can approximate the integral as 

_ j_j_ dn _ cw _ 

N ^ Jo dx cos^ 9 + dy sin^ 9 ^ 

VTE log ajv /(^)1 7 

7r3/2 jy TV 


OiN 


dx cos^ 9 + dy sin^ 9 
2y^cos2 9 + ^ sin^ 9 


+ 1 


where /(A) is a nonsingular function that depends on the cutoff and I = log (A/if) . 

Next we calculate Here in the numerator compensate the singularity of d{il,x,y), and hence the 

large-A^ limit can be safely taken before the integration. Then we obtain 




1 dfl f , , x'^ + y"^ dxx'^ 

Att^N fl J ^ ^{1 + x'^ + y"^)^ _(l + cx^ + y'^y/^ 


dyy"^ 

(1 + cx'^ -I- ?/2)3/4 


0.5561, 

N 


(46) 


The integral with respect to x and y is obtained numerically. 


2 . T,k 


We assume |n|, Eq with uj = ky = 0 and expand Go(H, -I- qx,qy) up to Using the dimensionless 

parameters x and y, we have 


'Eh = 


( 2 . 


with 




f dn f (1-bj/2)2 - 12x4(1-by2)-f 3a;® ^ _ ^.(o) , ^(i) 

j Wlj - (l+x7 + i,»P -^ + U’. 

3x® — 13x4 — 12x4y2 + 


(2^)3 y \n\ 


St / dxdy- 


(1 -I- x4 -I- 2/2)3 


-d{n,x,y). 


In calculating the large-AT limit can be safely taken similarly as in and we obtain 


Ai) 


1 


47r3AT 


ffdxdy 


3x® — 13x4 — 12x42/2 _|_ y4 _|_ ^2 
(1 -I- x4 -I- 2/2)3 


dyy"^ 


(1 -I-Cx4 -I-2/2)3/4 (1 -I-cx4 -I-2/2)3/4j 


(47) 


(48) 


0.4521. 


N 


(49) 


3. Eky 


Finally we consider the case where |U|, Eg ^ \vky\ with to = kx = 0 and expand Go(U, qx,ky + qy) up to the linear 
order in ky. Then we obtain 


Eh... = 


a 


(2^)3 7 |U| 


St / dxdy 


1 -f x4 — 2/^ 

(1 -I- x4 -I- 2/2)2 


d(U,x,2/) = S(0)+44) 


where the function is defined by 


= 




(2^)3 J |0| 

and calculated similarly as and as 


T7TT / dxdy 


4 2 

X — y 


(1 -I- x4 -I- 2/2)2 


1 


fcy 47r3AT 


r dn 

Ie ^ 


dxdy 


4 2 


(1 -b x4 - I - 2 / 2)2 


dxX^ 


d{n,x,y), 


dyy'^ 


(1 -b cx4 -b 2/2)4/4 (1 -(-c 2;4 _|_ y2)3/4 


(50) 


(51) 


0.2157, 

N 


(52) 
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III. RENORMALIZATION GROUP ANALYSIS 


Now let us perform a renormalization group analysis. The effective action is given by 

f '>pi{-'ii^ + AklTx + vkyTy-J:{k))tpk+ [ * 5 ( 1 + ^ [ D~^{k)(l)l(j)k 

Jk Jk,k' ^ Jk 

= *'^(1 + + ^(1 + '^kjklr^ + v{l + T,kJkyTy)ipk 

Jk 

+ [ ig{l + ST{k,k'))tjjlilJk'4>k-k'+ 1: f D~'^ {k)4)\(j)k. (53) 

Jk,k' ^ Jk 

The action changes by rescaling as 

5 f + s„) + Z^^b-^A{1 + ^kSklr. + + SfcJfc^rJ 

J k ^ ^ 

+ f [ D-\k)^l^-^, (54) 

J k,k' ^ J k 

where we put tildes to represent rescaled quantities. Here we rescale the frequency and momenta as 

uj = b^oj^ kx = bkx, ky = b^^ky, (55) 

and the fields and the other parameters are rescaled as 

■0 = 0 = ^00, A = ZaA, V = ZyV, g = Zgg. (56) 

We note that D~^ has the same scaling dimension as kx- The effective action keeps its form by the rescaling to be 

^ + + ^hTv)^k + ^ W^li’k'h-k' (57) 


Comparing Eqs. (54) and (57), we obtain the following relations: 


= 6-1-2"-^= (1 + E^), 
^2 = 6-2— 


ZA = b 

z„ = b 


— A-3-2-22 C7-2 


^l/> (1 "I” ~ 5^ (1 + ~ ^w) ~ 1 + (z — 2 + 7 a)0 


— 1 —2 —220 '7 — ^ 2 , 


^ib (l + ^fc„)~5^ '^^(1 + S/j — E^j) R:! 1 + (z — Z 2 + 7i,)0 


Zy = b 


Z - Z 2 . 


■2-2"-2"^202^-1 (1 + 5r(o,0)) «6("-"=)/2(i + 5r(o,o)-s^) « i + ^^z 


where we define 7 ^ and 7 „ as 


7A (Sfc^ 


_ 0.1261 _d 
1 = 0 ^ N ^ 


^ 0.3625 
/=o N ■ 


(58) 


(59) 


IV. CHARGE SCREENING 


Let us consider the static screening effect by introducing a charged impurity. The charge distribution induced by 
the impurity is given by 


/Oind(q) = ZeD{uj = 0, g)lVn(a; = 0, q) 

Bx 1 I + By \/|gy| 


= —ZeotA! 


where the static polarization is 


n(" = “'<J) = -i6l«-l-60fr(3/4) 


27</^ + </y + (Sylg.-I + By^/\qbf 

a r(5/4) 


(60) 


^kyl ^ — Ot(^x\qx\ CX.By^ 1% 1^ • (61) 
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FIG. 5: Induced charge distribution pind(r’). x and y are measured in unit of {B^/By)^ and we set ZeNaB^ ^Byf2 = 1. 


The real space distribution of the screening charge /9ind(f) can be obtained after Fourier transformation: 

Pind(r) = J ^^^e*'» ''pi„d(q). (62) 

We define partially integrated charge densities as 

Qxix) = J dypindir) = J = 0 )> 

Qyiv) = j dxpindir) = j ^e^’^^^Pindiqx =0,qy)- (63) 


A. ajv <C 1 

For A^a <C 1, the induced charge in momentum space is given by 

Bx 19a: I + By \/|(7y I 


Pindiq) = -ZeNa- 


and hence we obtain the partially integrated charge density 


‘2Jql + qi 


Qx{x) = -l-ZeNaBxS{x), Qy{y) = - ^ZeNaBy—^. 

2 VStt ^\y\ 


(64) 


(65) 


Thus the integrated charge along the quadratic dispersion direction Qx{x) is localized at the impurity site whereas 
the integrated charge along the linear dispersion direction Qy{y) has a power law decay. 

The charge distribution in real space /9ind(’’) is obtained by the Fourier transform of Pind(9): 


And(g) = ^ZeNa[Bxhi{r) + Byh2{r)], 


where hi[r) and h 2 {r) are defined as 


h 2 {r) = - 


hi{r) = 


1 


1 - 3^sinh"^ A 


'(h) 


v/27r3/2|a;|3/2 l + p2 


Re 


2E 


TT^ r 

l+ip 
2 


- K 


l + ip 
2 


— pim 


K 


l+ip 

2 


( 66 ) 


(67) 


where K{m) and E(m) are the complete elliptic integrals of the first and second kind, respectively, with the parameter 
m, and we define r = yjx'^ + y"^ and p = |j//x|. The charge distribution is shown in Fig.j^ 
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B. ajv ^ 1 

For Na » 1, Pindig) = —Ze leads to Pind(<7) = —ZeS{r)l{2Tr)‘^. The effect of the screening is large in this limit, 
and induced charge is localized on the impurity. 


V. PHYSICAL OBSERVABLES 

If we consider an observable O, we expect the scaling behavior 

0{k,uj,v,A) = zoO{k{l),uj{l),v(l),A{l)), (68) 

where zq is the scaling dimension of the operator O, and k{l) = {zx{l)kx,e^ky) and = Zuj{l)oj are the running 
momenta and frequency, respectively. Zx and z^j obey 

Defining b = e\ we obtain 

Zx = z^ = h^-^-. (70) 

The scaling behavior discussed here is only valid for ^ 1. In the weak coupling case, the low-energy physics can 
be obtained via naive Hartree-Fock perturbation theory with weak coupling RG analysis. Here one expects logarithmic 
corrections, not corrections to power law behavior. 


A. Compressibility 


In case of the particle density O = n, we have zo=n = b such that the compressibility k = dn/dp follows 

zo=K = b~^Zx^Zui. The scaling relation then implies 


k{T) = 6-T'''-(i+T''^-T'’')/2K(6i-'y“r), 


(71) 


which leads to 


k(T) CX rl/2+‘^, 


with 


</> = 7 « + 


0.4255 

N 


(72) 


(73) 


B. Heat capacity 


The scaling of the free energy density is 

/(T, p) = b-^Zx^zZ^f{z^T, z^p), (74) 

which reproduces the above scaling dimensions for the particle density via n = df/dp. For the heat capacity 
C = df /dT, its scaling dimension implies 


C{T)=b-^Zx^C{z^T). 


(75) 


Then we obtain 


G(r) oc t3/2+<^. 


(76) 
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C. Optical conductivity 


We first calculate first the bare optical conductivity. The current-current correlation function liafiiq, w) {a, P = x,y) 
for the two direction is 


Ilxx{i^) = -ANA^e^ J ^^^^kltT[TxG{u}, k)TxG{u: A O, fc)], 

Ayy{i^l) = -Nv^e^ f ^ tT[TyG{u;, k)TyG{uj -I- O, k)]. 


(77) 

(78) 


Since we are only interested in the imaginary part of IIqq on the real axis, and thus we subtract the 0 = 0 contribution. 
Then we obtain 


n..(*o) n,.(o)- 

Uyyiin) - nj^y(o) = -Nn^v^e^ J 


d^k 


Akl • vX 


1 


( 2^)2 ( yl 2 fc 4 + „ 2 ^ 2 ) 3/2 ^2 + 4 ^ 2 fc 4 + 4 ^ 2^2 ’ 

dPk A^kt _1_ 

(27r)2 {A'^k'^ + v‘^ky)^G 02 -|- 4A‘^k^ + Av'^ky 
With X = VAkx/^/l^l and y = vky/\^\, it follows that 


(79) 

(80) 


Uxxm - n,,(0) = -4c,iV|0|3/2^e2, 

V 

n,,(zO) - n,,(o) = -cyN\n\^/^^e^ 


where and Cy are numerical constants, given by 


Cx — 


dxdy 


2 2 
X y 


1 


(27r)2 (a;4 -|- y2^3/2 2 ^ _|_ 42;4 _|_ 4 y 2 ’ 


dxdy x'^ 1 

J_ao (27r)2 (x^ + 1 -|- 4x4 -|- 4j/2 


(81) 

(82) 


(83) 


It holds for —1 < ?7 < 1 that nQQ(zO) = —i3|0|^ leads to ImnQQ(w -|-10+) = B cos(7r?7/2)a;|w| Thus, it follows 

for the real part of the optical conductivity 


that 


fJcKQ, (^) 


ImnQQ,(a; + f04") 

UJ 


^ XX 


(w) =fV 




ayy{uj) =N—c' 



(84) 


(85) 


with Wo = /A and the constants = AcxI'J^ ss 0.05, c’y = ~ 0.05. 

Using the gauge invariance, one hnds that the scaling dimensions of the optical conductivity are different for the 
two directions: 


Then it follows that 


crxx{‘^) = b Zx(Jxx{z^^), CTyy{uj) =bz^^ayy{z^uj) 


axxii^) oc N 


g2 / ^ X 1/2+4:. 


h Vwo 


ayy{uj) oc N 


^2 . X -l/2-0„ 


h Vwo 


( 86 ) 


(87) 


with (j>cr = Jv — 7a/ 2 « 0.299/A^. The anisotropy of the optical conductivity is amplihed by strong interactions. 


D. Diamagnetic susceptibility 


We consider the diamagnetic response to a field perpendicular to the plane. We use for the diamagnetic susceptibility 


Xd = lim- Ilxyiq,uj = 0), 

q^O QxQy 


(88) 



15 


with the current-current correlation function Ilai 3 {q,uj) {a,j3 = x,y). The scaling dimension of IVxy is z~^. This 
follows from the following logic: Using the Kubo formula and our results for the conductivity it holds that Ilxx has a 
scaling dimension b~^ZxZ~^ and IVyy has . Thus, a mixed term lixy must have the geometric mean as scaling 

dimension, i.e., z~^. It leads to the scaling relation 

Xd{T) = hZxZ~'^XD{zu,T), (89) 


which yields 


XD oc 


rjn—1/2 —(p 


(90) 


with (f> = Xv + IaI‘2' = 0.4255/-/V. 



